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Abstract Integrable operators arise in random matrix theory, where they describe the 

asymptotic eigenvalue distribution of large self-adjoint random matrices from the general- 

^ | ized unitary ensembles. This paper considers discrete Tracy-Widom operators, and gives 

sufficient conditions for a discrete integrable operator to be the square of a Hankel matrix. 

Examples include the discrete Bessel kernel and kernels arising from the almost Mathieu 

^ ■ equation and the Fourier transform of Mathieu's equation, 
-i— > 
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1. Introduction 

We consider Tracy-Widom operators arising from first-order recurrence relations 



a(j + l)=T(j)a(j) (j = l,2,...), (1.1) 



ON 

where a(j) is a real 2x1 vector and T(x) is a 2 x 2 real matrix with entries that are 
rational functions of x, and such that detT(j) = 1. Then with the skew-symmetric matrix 



J 



-1 

1 



the Tracy-Widom operator K : i (N) — > I (N) has matrix 

K{m,n) = (m 7^ 77,) (1.2) 

m — n 

with respect to the usual orthonormal basis. In specific examples, there are natural ways 
of defining the diagonal K(n,n), as we discuss below. We recall that a real sequence 
((f) j) G £ 2 (N) gives a Hankel matrix = [<f>j+k-i]j°k=v c l earr y ^<f> is symmetric, and 
is Hilbert-Schmidt if and only if ^2^ = iTi\(f) n \ 2 converges. See [14]. 

We consider whether such a K may be expressed as K = T^, where is a Hankel 
matrix that is self-adjoint and Hilbert-Schmidt. 
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A significant example from random matrix theory is the discrete Bessel kernel 



B(m n) = ^ Jm(2v^)J n -i(2y^) - J n (2y^)J m _ 1 (2^) 



m — n 



(1.3) 



as considered by Borodin et al [5] and Johansson [11]. They showed that B is the square of 
the Hilbert-Schmidt Hankel matrix [J m+ / C _i(2v / ^)], and thus obtained information about 
the spectrum of B itself. 

Tracy and Widom observed that many of the fundamentally important kernels in 
random matrix theory have the form 



W(x,y)= m9{y) - fiy)9{x) (x^y) 



x-y 



(1.4) 



where /, g are bounded real functions in L 2 (0, oo) such that 

d_ 

dx 



7W" 




a(x) 


(3{x) ' 




"/(*)" 


_g(x) _ 






—a(x) 




g{x) 



(1.5) 



and a(x), (3{x) and 7(2;) are real rational functions; see [16]. Then there is a bounded linear 
operator W : L 2 (0, 00) — > L 2 (0, 00) given by Wh(x) = W{x,y)h(y)dy. Of particular 
importance is the case in which W is a trace-class operator such that < W < I, since 
such a W is associated with a determinantal point process as in [15]. In order to verify 
this property in special cases, Tracy and Widom showed that W is the square of a Hankel 
operator V which is self-adjoint and Hilbert-Schmidt. The spectral theory and realization 
of Hankel operators is well understood [8, 13], so such a factorization is valuable. In [3, 
4], we developed this method further and considered which differential equations lead to 
kernels that can be factored as squares of Hankel operators. The factorization theorems 
involve the notion of operator monotone functions. 

Here we consider various discrete Tracy-Widom kernels and their factorization as 
Hankel products, using the formal analogy between differential equations and difference 
equations which suggests likely factorization theorems. We can write the matrix in (1.5) 
as JO(x), where Q(x) is real and symmetric, and then consider the analogous one-step 
transition matrix to be T(x) — exp( JQ(x)). The functions x ^ x and 1 — >■ — 1/x are operator 
monotone monotone increasing on (0, 00), and they appear in the transition matrices for 
the discrete Bessel kernel in section 3 and the discrete analogue of the Laguerre kernel in 
section 4. 

The functions x ^ x 2 and x 1— > — 1/x 2 are not operator monotone increasing by 
[9], and so we cannot hope to have simple factorization theorems when they appear in 
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the transition matrix T. In the special case of the parabolic cylinder equation — 4>'^{x) + 
(x 2 /4 — l/2)(f) n (x) = n(f) n (x) : Aubrun [1] recovered a factorization of the corresponding 
kernel 

K , x x = <Pn{x)<Pn-l(y) ~ (pn-l(x)<f> n {y) 

x-y 

in the form 

K = T^T i ,+T i ,T^ (1.6) 

where and are bounded and self-adjoint Hankel operators. The parabolic cylinder 
function is a confluent form of Mathieu's functions associated with the elliptic cylinder, 
since the parabola is the limiting case of an ellipse as the eccentricity increases to one [17, 
p. 427]. Hence it is natural to factorize kernels associated with Mathieu's equation 

d 2 u 

— -j + (a + (3 cos z)u(z) = (1.7) 

in the form (1.6). In section 5 we consider the first-order difference equation associated 
with the Fourier transform of Mathieu's equation. In section 6, we consider almost Mathieu 
operators. 

For a compact and self-adjoint operator W, the spectrum consists of real eigenvalues 
Xj which may be ordered so that the sequence of singular numbers Sj = \Xj\ satisfies 
si > s 2 > — While the factorization W = V 2 immediately determines the spectrum of 
W from the spectrum of T, a factorization (1.6) imposes bounds upon the singular numbers 
of K is terms of the eigenvalues of and T^. 

2. Factoring discrete Tracy— Widom operators as squares of Hankel matrices 

Theorem 2.1. Let T(j) and B(j) be 2 x 2 real matrices, and (a(j)) a sequence of 
real 2x1 vectors such that 

a(j + l)=T(j)a(j) (j G N), (2.1) 
a(j)-0 (j-oo), (2.2) 

oo 

J2\\B(j)a(j)\\ 2 <oc. (2.3) 

Suppose further that there exists a real symmetric matrix C with eigenvalues and X, 
where X < 0, such that 



T{nfJT{m) - J 



m — n 



= B(nfCB(m) (m ^ n; m, n G N). (2.4) 
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Let 4>{j) = \\\ l / 2 (v\,B{j)a{j)) where v\ is a real unit eigenvector corresponding to 
A. Then is compact and K = T| has entries 

(Ja(m),a(n)) , . . 
K(m 7 n) = ± — v h \ )i m ^n;m,nGN. 2.5 
m — n 

Proof. Let K(m 7 n) be as in (2.5) and let 

oo 

G(m, n) = K{m, n) - y~] <j>(m + k- l)(j>(n + k - 1), (m ^ n; m, n G N) (2.6) 

fe=i 

where the infinite sum converges because of the condition (2.3). Observe that a(j) — > 
implies that the first term in G(m, n) tends to zero as m or n — > oo, and that the same is 
true of the Hankel sum: 

oo / oo \ V2 / oo X 1/2 

|5^(m + fc-l)0(n + fc-l)| < ^|(/)(m + A;-l)| 2 l</>(™ + ^ " 1) N -> 0, 

fc=i \fe=i / \fc=i / 

(2.7) 

as m or n — > oo. Hence G(m, n) — > as m or n — > oo. Now let U be the real orthogonal 
matrix with v\ in the first column, and the eigenvector of C corresponding to the eigenvalue 
in the second column. Then U t CU = diag(A, 0), and we have 

K(m + 1, n + 1) — K(m, n) = — (JT(m)a(m),T(n)a(n)) - (Ja(m),a(n)) 

! ((T(nyjT(m) - J)a(m),a(n)) 



m — n 

= {B^nfCB^a^m^a^n)) 

= (C/diag(A, 0)*7'B(m)a(m), B(n)a(n)) 

= A <diag(l, 0)t/ t S(m)a(m), diag(l, 0)f7*.B(ra)a(ra)) 

= -(j>(m)(j>(n). (2.8) 

The above calculation, and the equality 

oo 

^2 <P( m + k)(f)(n + k) - <p(m + k- l)<j>(n + k - 1) = -(j>(m)(j>(n) (2.9) 
k=i 

together imply that G(m+1, n + 1) = G(m, n), and so in fact G(m, n) = for all m, n G N, 
which gives 

oo 

K(m, n) = (j)(m + k - l)(j)(n + k - 1) (m 7^ n; m, n G N), (2.10) 
fc=i 



where the right-hand side is the (m, n) th entry of the square of a Hankel matrix Y^. Further, 
we observe that K is the composition of the Hilbert transform with the compact operators 
^ 2 (N; C) — > £ 2 (N; C 2 ) given by (x n ) h- > (Ja(n)x n ) and the adjoint of (x n ) h- > (a(n);r n ); so 
if is compact. We deduce that is also compact. 

□ 

Proposition 2.2. Let a(j), T(j) and B(j) satisfy conditions (2.1) and (2.2), and 
suppose further that 

OO 

J^jWBUHM 2 <oo. (2.11) 

Now let K = be as in Theorem 2.1. 

(i) Then K is a positive semideSnite and trace-class operator. 

(ii) For each n G N, there exist self-adjoint Hankel operators Y n , where Y n has rank 
at most n, such that 

Sn (K) = ||rv-r n || 2 , (2.12) 

so T 2 — > K as n — > oo. 

Proof, (i) The Hilbert-Schmidt norm of Y^ satisfies HT^H 2 ^ = YlkLi k(/)(k) 2 < oo. 
Hence K = T 2 is of trace class. 

(ii) Since Y^ is self-adjoint, the singular numbers satisfy s n (K) = s n (T 2 ) = s n (r^) 2 . 
By the Adamyan-Arov-Krein theorem [14] , there exists a unique Hankel operator F n with 
rank at most n such that s n (Y ( j } ) = ||r^ — r n ||. Evidently T* is also a Hankel operator of 
rank at most n such that s n (r^,) = ||r^ — T* ||, so by uniqueness Y n = T*. 

We have ||r n — r^,|| — >• as n — > oo, so r 2 — * T 2 as n — > oo. 

□ 

Definition. For a compact and self-adjoint operator If on a Hilbert space H, the 
spectral multiplicity function vyy : R — >• {0, 1, . . .} U {oo} is given by 

v w (\) =dim{x e H :Wx = \x} (A G R). (2.13) 

Proposition 2.3. Let K be as in Proposition 2.2. Then the following hold: 

(i) u K (0) = or u K (0) = oo; 

(ii) z/^(A) < oo and v K (A) = ^(v^) + vr^—y/X) for all A > 0. 
(Hi) If vk{X) is even, then vtAVX) = vr+i—VX)- 

(iv) If v K (X) is odd, then \vr^(y/X) - ^(-\/A)| = 1. 

Proof, (i) follows from Beurling's theorem (see [14, p. 15]), while (ii) is elementary. 
Peller, Megretskh and Treil show in [13] that for any compact and self-adjoint Hankel 
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operator T^, the spectral multiplicity function satisfies |^r^(A) — z/ r<#) (— A)| < f . Using this, 
and (ii), statements (Hi) and (iv) follow immediately. 

3. The Discrete Bessel Kernel 

We show how Theorem 2.1 can be applied to the discrete Bessel kernel to recover a result 
from [5, 11]. 

Proposition 3.1. Let J n (z) be the Bessel function of the first kind of order n, let 
Jn = Jn(2\/0), where 9 > 0; let 4>(n) = J n +i and a(n) = [v^Jn, Jn+i]*. Then the Hankel 
operator is Hilbert-Schmidt, and B = V 2 , has entries 

B(m,n)= ( Ja (m),a(n)) ^ m _^ n . m ^ neN y (31) 
m — n 

Proof. It is clear that (2.1) holds, since we have the recurrence relation 

n 4- 1 

J n+2 (2z) = J n+ i(2z) - J n (2z), (3.2) 

z 

giving a(n + 1) = T(n)a(n), where 

T(n) = 

Note that 



o Ve 

-1 n+l 



(3.3) 



T(nYJT {m) - J = 
m — n 

where C = diag(0, — 1), which is clearly of rank one. The non-zero eigenvalue of C is 
A = —1, and a corresponding unit eigenvector is v\ = [0, 1]*, so 

|A| 1/2 (fA,a(n)) = J n+ i = 4>(n). 

We now verify condition (2.11), and thus (2.2). Note that 

1 1 x - 

n=l n=l n=l n=l 

The standard formula (see [17, p. 379]) 

oo oo 

e i2Vflsin^ = j (2v^) + 2 J2m(2v^)cos2m^ + 2z J2m-i(2v^)sin(2m- 1)^ (3.6) 



m=l m=l 



and Parseval's identity can be used to show that J (2v^) 2 + 2 YZ=i Jm(2v^) 2 = 1 for all 
9 > 0, and hence that the sum on the right hand side of (3.5) is finite. 

4. A discrete analogue of the Laguerre differential equation 

In this section we consider a case in which condition (2.2) is violated and we cannot hope 
to factor the kernel K as the square of a Hankel operator. Nevertheless, we can identify a 
Toeplitz operator W such that K — W factors as a product of Hankels. 

Proposition 4.1. For 9 G R, let (a(j)) satisfy the recurrence relation a(j + 1) = 
T(j)a{j) with 

~ 0/(3 + 1) -1' 



1 



(4.1) 



and a(l) = [6>, 1]*. Then there exist polynomials Pj(9) of degree j such that 

(i) a(j) = [p J (9),p J . 1 (9)] t . 

(ii) The self-adjoint Hankel matrix = [<j)(j + k — l)]°° fc=1 with entries 

m = f& (4.2) 

is a bounded linear operator such that 9T^ = K + W where K has entries 

tw \ (Ja(m),a(n)) . . _ TS oN 

K(m,n) = - — v ' v " (m ^ n;m,n G N), (4.3) 
m — n 

and W is a bounded Toeplitz operator with matrix 

Vr (m — n) = (m n ) 

m — n 

for some 2x2 matrix T^, and m and n are the congruence classes of m and n modulo 4. 

Proof. When 9 = 0, the recurrence relation reduces to a(n+l) = Ja(n), with solution 
a(n) = J n ~ 1 a(l). This gives rise to a kernel 

= (Jn ( ro ), n(n )) = (J "^-.(l), a( l)) 
m — n m — n 

which has the shape of a Toeplitz operator, and is a variant on the Hilbert transform 
matrix [1/ (m - n)] m ^ n . 

Now suppose 9^0. The matrix J satisfies J 4 = /, and so we consider the partial 
product of the T(j) in bunches of four, with the j bunch giving 

B(j) = T(Aj)T(Aj - l)T(4j - 2)T(4j - 3) = J - A J + 0(l/j 2 ) (j G N). (4.5) 

2j 



Now we deduce that 

\\B(j)f = \\B(jyB(j)\\ = 1 + 0(1/ f), (4.6) 
and likewise with B(j)~ x in place of B(j); so there exists C(9) such that 

\\T(n)T(n-l)...T(2)T(l)\\<C(9), 

||T(l)- 1 T(2)- 1 ...T(n-l)- 1 T(n)- 1 || < C(9) (n = l,2,...). (4.7) 

It follows that there exists k(9) > such that k(9) < ||a(n)|| < k(9)~ x for all n, so (2.2) is 
violated. We introduce 

k 1 

C k = e W (ejJ2^)B(k)B(k - 1) . . .B(l), (4.8) 
i=i 3 

which satisfies C^+i — Ck = 0(1/ k 2 ); so the limit 

= lim C k (4.9) 

fc^oo 



exists. One can check that 

k+n „ T k+m Q , 



i=i J i=i J 

as /c — > oo, and hence 

(Ja(m + 4/c),a(n + 4k)) 

= (JT(m + 4k)T(m + 4k - 1) . . . T(l)a(l), T(n + 4k)T(n + 4k - 1) . ..T(l)a(l)} 
-> (r^T^afl), J^TooaCl)) (A; -> oo). (4.11) 

For temporary convenience we introduce 

K(m, n) = \ (m-n) tor m ^ n ' (4.12) 
I for m = n. 

The discrete Hilbert transform is bounded on £ 2 (Z;C 2 ) by [7], so K defines a bounded 
linear operator on ^ 2 (N), but condition (2.2) is violated. 
The pj (9) satisfy the recurrence relation 

Pn+1 (9) + Pn - 1 (9) = -?—Pn(9) (4.13) 

n + 1 



with po(9) = 1 and pi(9) = 0, so clearly Pj(9) is a polynomial of degree j such that 
\Pn{9)\ < k(^)" 1 for all 9 and n. Further, 



T(nyjT(m) - J 



m — n 



-0/(m+l)(ra+l) 




(4.14) 



so that, by the calculation in the proof of Theorem 2.1, 



K(m + 1, n + 1) - K(m, n) = -9 



Pm{9)Pn{9) 

(m + l)(n + l) 



(4.15) 



We can write 



TV 



K(m, n) - K(m + N + l,n + N + 1) = ^ 



fc=0 



Pm+k{9)p n +k{9) 

(m + k + l)(n + k + 1) 



(m^n) (4.16) 



where the limit 



lim K(m + N + 1, n + N + 1) = W(m - n) (m 7^ n) 
exists and is finite since the sequence 4>(j) — Pj(9)/(j + 1) is square summable; so 



(4.17) 



K(m, n) = ^ — Pm + k ^ Pn ^\ ; : + ir(/// - (/„ / „). ( 1.18) 



, o (m + k + l)(n + k + l) 
We now define W(0) =0, and let 

TS( \ flV^ Pm+k(9)Pn+k(9) , . tA-\C\\ 

A m,n =0> - — - 1 — - + W(m — n): 4.19 

f4 {m + k + 1 ){n + k+l) y h y ! 

k=0 

so that the matrix of K equals the matrix of K, except on the principal diagonal, and the 
principal diagonal of K is a bounded sequence; hence K is a bounded linear operator and 
satisfies the preceding identities also for m = n. Let £ be the shift operator on £ 2 (N). 
Now 9T^ equals the limit in the weak operator topology of the sequence K — S* n KS n as 
n — > 00, so is bounded and hence W = K — is also bounded. We recognise the 
matrix of W as 

un , (J^+ 1 T 00 q(l),T 00 q(l)) 

W{m — n) = (mf^n). 

m — n 



□ 
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Remark 4.2. The generating function f(z) = J2jLoPj(@) zj satisfies the differential 
equation 

(l + z 2 )f'(z) + (2z-6)f(z)=0 



with initial condition /(0) = 1, and hence 

l-iz\M/ 2 1 



l + iz 



1 + z 



2 ' 



(4.20) 



For comparison, Laguerre's equation [16] may be expressed as 



d 


u(x) 




r 




u(x) 


dx 


u'{x) 




l/4-(n+l)/x 




u'(x) 



(4.21) 



with solution u(x) = xe x / 2 Ln (x) where 



x 1 e x d n 
n\ dx n 



(x n+1 e- x ) (x > 0) 



(4.22) 



is the Laguerre polynomial of degree n and parameter 1. The Laplace transform of u is 
the rational function 



C{u; A) = (n+ 1) 



(A 



(A+|)-+2 



(»A > -1/2). 



(4.23) 



5. The Fourier transform of Mathieu's equation 

Let P 7^ be a real number; then there exists a sequence of real values of a such that 
Mathieu's equation 



d^u 

— + (a + /3 cos e)u(9) = 
dt) z 



(5.1) 



has a real periodic solution with period 2ir or 4rr. The odd or even periodic solutions 
are known as Mathieu functions, and various determinants describe the dependence of the 
eigenvalues a on /?, as in [12, 17]. Here we are concerned with some matrices that arise 
from the Fourier transform of the differential equation. 

Theorem 5.1. Suppose that u has Fourier expansion u{9) = Xl^L-oo b n e tne . Let T u 
be the Hankel matrix [bj+k-i}fk=i' ^ be the Hankel matrix [(j + k — l)bj + k-i]j°k=i 
and let K = (—2/ f3)(T u T v + T V T U ). Then K is a trace class operator on £ 2 (N) such that 



trace(-ftT) 



1 



du i 



(3% J l d6 



d9 



and 



K(j, k) = b Jz2hzMh=l 
j — k 



(5.2) 
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Proof. Since u is real, we have b m = b- m . The recurrence relation for the Fourier 
coefficients 

2(-n 2 + a)b n + [3b n+1 + [3b n -! = (5.3) 



may be expressed as the first-order recurrence relation 



f 

-f (2//3)(re 2 -a) 



b n 




_b n +i _ 





b n -l 
b n 



(5.4) 



or in the obvious shorthand a (re + 1) = T(re)a(re). Then we have 

T{nYJT{m) - J 



m — re 





(-2//3)(m + re) 



(5.5) 



We introduce the kernel K by the formula 

( Ja(m), a(n)) 



bnbm-i — b n -ib r 



m — n 



m — re 



K(m 7 re) 
which therefore satisfies 

K(m + l,re + 1) - K(m,n) = (-2/ (3)(mb m b n + nb n b m ) 
and K(m, re) — > as m, n — > oo in any way such that m 7^ re. We deduce that 
2 ^ 

K(m, re) = - 2_^{m + k)b m+k b n+k + (re + k)b n+k b m+k (m, re G N; m 7^ re). 



(5.6) 



fc=0 



Hence K(m, re) is the (m, re) entry of the matrix of K = (2/ f3)(T u T v J r Y v Y u ) for all m 7^ re. 

Since u and it" are square integrable, the series Xl^L-oo n4 |^| 2 converges; so T u and 
r„ are Hilbert-Schmidt, and K is trace class. Further, we have 

00 

trace(-ftT) = ^ K(m, m) 



m=l 



m+k — 1 



m,k=l 



00 



2^2 
m 

m=l 

2 <tf 



(5.7) 
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6. Almost Mathieu operators 

We introduce the almost Mathieu operator H : £ 2 (Z) — > £ 2 (Z) by 

{Hu) n = u n+ i + u n -i + A cos 2n(n9 + a) u n (6.1) 

for u = (u n )'£L_ 00 G ^ 2 (Z), where (Hu) n denotes the n th term in the sequence Hu G £ 2 (Z). 
For all real A, 9, uj and a, the operator H is bounded and self-adjoint, with spectrum 
contained in [—2 — |A|,2 + |A|]. According to the precise values of the parameters, as 
we discuss below, the spectrum can consist of a mixture of point spectrum, continuous 
spectrum and singular continuous spectrum. 

Definition. Let E be an eigenvalue of H with the corresponding eigenvector (u n ). 
Say that (u n ) has exponential decay if there exist C, 5 > and no G Z such that 

\u n \ < C e - 5 l n - n °l (n G Z). (6.2) 

(Typically, 5 depends upon 6, a and A in a complicated fashion.) Say that H exhibits 
Anderson localization if its spectrum is pure point and all eigenvectors are of exponential 
decay. 

Definition. Say that 9 G R is Diophantine if there exist c{9) > and r(9) > such 

that 

|sin27r^|> C (^)|j|- r(e) C?GZ\{0}). (6.3) 

With respect to Lebesgue measure, almost all real numbers are Diophantine; see [2, 
p. 373]. Clearly rational numbers are not Diophantine, nor are Liouville numbers; see [10]. 

Following the work of several mathematicians, as summarized in [6], Jitomirskaya [10] 
obtained a satisfactory description of the spectrum of the almost Mathieu operator. 

Lemma 6.1. [10] Suppose that 9 is Diophantine. Then there exists a set Sq such 
that R \ Se has Lebesgue measure zero, and such that for all A > 2 and all a G Se, the 
Mathieu operator H has an eigenvalue E such that the corresponding eigenvector (u n ) is 
of exponential decay as n — * ±oo. 

Moreover, she proved that the Mathieu operator has pure point spectrum for A > 2, 
and further conjectured that the same conclusion holds for all real a. For comparison, for 
A = 2 there is purely singular continuous spectrum; whereas for < A < 2 there is purely 
absolutely continuous spectrum. 

We introduce the Hankel matrices 

r c =[cosir(a + x0 + k6)u x+k ]% k=1 , T s = [sin n (a + x9 + k9) u x+k ]^ k=1 . (6.4) 
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There exists a bounded and measurable function $ : T — > M 2 (C) such that 



cos 7r(a + kO) Uk smn(a + k9) Uk 
cos 7r(o! + k9) Uk cos 7r(ct + k6) Uk 



(6.5) 



the block Hankel operator associated with $ is [<l(j + fc)], which becomes, after a rear- 
rangement of the block form, the matrix 



r c r s 



(6.6) 



The negative Fourier coefficients of $ are not uniquely determined by but may be 
chosen advantageously. Note that all of these operators are self-adjoint. 

We introduce operators K and L by 



L 


K 






-T(Ts ~\~ r s r c 


K 


L 









= r 



(6.7) 



Theorem 6.2. Let u n be as in Lemma 6.1. 
(i) The matrices of K and L are given by: 



K(m, n) = 



U m —\U n U n — \U r 



(m 7^ n; m, n e N); 



2A sin7r6>(m — n) 

oo 

L(m, n ) = cos 7r9(m — n) u m +kUn+k (tti, n G N). 



(6.8) 



fc=i 



(j'i) Then T| is a positive semidehnite trace-class operator. 

(in) Further, the eigenvectors of K, T c and T s are of exponential decay as x — > oo. 
(iv) There exists a bounded and measurable function \I> n : T — > M 2 (C) such that the 
associated Hankel operator [& n (j + k)] has rank less than or equal to n and 



Sn(r$) 



-r 



1$ - ^ T 



(6.9) 



(v) The eigenvalues of K, L, T s and T c are of exponential decay. 

Proof, (i) First we observe that the formula for K(m, n) makes sense for m ^ n since 
9 is irrational. The discrete Mathieu equation 



u n +i + u n -i + A cos 2n(n9 + a)u n = Eu r 



(6.10) 
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gives the system 







_U n+ l _ 





1 

-1 E — Xcos2n(n6 + a) 



Writing T(n) for the one-step transition matrix, we have 
T(n)VT(m) - J 



U n -1 

U„. 



(6.11) 




1 



2Asin ix9{m — n) 
so with a(m) = [it m _i,it m ]* we introduce 



(sin7r(n6' + a) cos7r(m6> + a) + cos n(n6 + a) sin7r(m6> + a)) ; 

(6.12) 



K(m, n) 



( Ja(m), a(n)) 
2A sin7r6>(m — n) 



(6.13) 



which satisfies 



K(m + 1, n + 1) — if (m, n) = (sin 7r(n6> + a) cos n(m6 + a) 

+ cos7r(n6> + a) sin7r(m6> + a))u n u m 

and if (m + /c, n + k) — > as — > oo. Hence by comparing entries of the products, we find 

that k = r c r s + r s r c . 

(ii) This is clear, since the entries of the matrix of T$ are real and summable. 

(iii) For any unit vector ip, we have 



T c <f(x) = ^2 cos n ( a + Ox + 9k) u x+ k^k 



(6.14) 



fc=0 

so by the Cauchy-Schwarz inequality we have the uniform bound 

1/2 

<) 

k=x 



00 1/2 

|r>(^)| < Q>> 2 fc ) , 



(6.15) 



where the right-hand side decays exponentially as x — > oo. A similar result applies with 
T s , so in particular the eigenvectors of T c and T s are of exponential decay at infinity. 

Now let ((pj)j%:i be an orthonormal basis of £ 2 (Z) consisting of eigenvectors of T s with 
corresponding eigenvalues aj. Then 



T c T s ip(x) = ^2a j {(p j ,(p)T c (p j (x) 



(6.16) 
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where 

oo OO -jy2 00 \j<2 

^2\<Tj{<Pj,<p)\ < ($^(v^i) 2 ) = HrjcHMI^- ( 6 - 17 ) 

This and a similar result for T s r c imply \Kip(x)\ decays exponentially as x — > 00. 

(iv) This is immediate from the vectorial form of the matrical AAK theorem [14] . 

(v) Since u n decays exponentially asn^ 00, we can approximate the Hankel matrices 
with finite matrices up to exponentially small error terms. For instance, we can approxi- 
mate T s by = [sin7r(a + x6 + kO)u x +kl{( x ,k):x+k<N}] which has rank less than N + 1 
and the operator norm satisfies 

00 

l|r s -rf)||< ( 6 - 18 ) 

k=N+l 

The s-numbers satisfy 

= s n (r|) = s n (r*) 2 . (6.19) 

□ 

A vectorial Hankel matrix has finite rank if and only if it has a rational symbol 
with coefficients of finite rank by [14, p. 19]. Peller provides a formula for the rank in 
terms of the coefficients. 
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